Kondo effect induced by a magnetic field 
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We study peculiarities of transport through a Coulomb blockade system tuned to the vicinity of 
the spin transition in its ground state. Such transitions can be induced in practice by application 
of a magnetic field. Tunneling of electrons between the dot and leads mixes the states belonging 
to the ground state manifold of the dot. Remarkably, both the orbital and spin degrees of freedom 
of the electrons are engaged in the mixing at the singlet-triplet transition point. We present a 
model which provides an adequate theoretical description of recent experiments with semiconductor 
quantum dots and carbon nanotubes. 
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I. INTRODUCTION 

Quantum dot devices provide a well-controlled object 
for studying quantum many-body physics. In many re- 
spects, such a device resembles an atom imbedded into 
a Fermi sea of itinerant electrons. These electrons are 
provided by the leads attached to the dot. The orbital 
mixing in the case of quantum dot corresponds to the 
electron tunneling through the junctions connecting the 
dot with leads. Voltage Vg applied to a gate - an elec- 
trode coupled to the dot capacitively - allows one to con- 
trol the number of electrons N on the dot. Almost at any 
gate voltage an electron must have a finite energy in or- 
der to overcome the on-dot Coulomb repulsion and tunnel 
into the dot. Therefore, the conductance of the device is 
suppressed at low temperatures (Coulomb blockade phe- 
nomenon iQ]). The exceptions are the points of charge 
degeneracy. At these points, two charge states of the 
dot have the same energy, and an electron can hop on 
and off the dot without paying an energy penalty. This 
results in a periodic peak structure in the dependence 
of the conductance G on V^. Away from the peaks, in 
the Coulomb blockade valleys, the charge fluctuations are 
negligible, and the number of electrons N is integer. 

Every time N is tuned to an odd integer, the dot must 
carry a half-integer spin. In the simplest case, the spin is 
S = 1/2, and is due to a single electron residing on the 
last occupied discrete level of the dot. Thus, the quan- 
tum dot behaves as 5* = 1/2 magnetic impurity imbedded 
into a tunneling barrier between two massive conductors. 
It is known |^ since mid-60's that the presence of such 
impurities leads to zero-bias anomalies in tunneling con- 
ductance [||, which are adequately explained Q in the 
context of the Kondo effect The advantage of the 
new experiments Q is in full control over the "magnetic 
impurity" responsible for the effect. For example, by 
varying the gate voltage, N can be changed. Kondo ef- 
fect results in the increased low temperature conductance 
only in the odd-A^ valleys. The even-A'^ valleys nominally 
correspond to the = spin state (non-magnetic impu- 



rity), and the conductance decreases with lowering the 
temperature. 

However, the enhancement of the low temperature con- 
ductance in the even-A^ valleys has been observed re- 
cently in both vertical and lateral [|j GaAs quantum 
dots, and in transport through carbon nanotubes Un- 
like the real atoms, the energy separation between the 
discrete states in a quantum dot is fairly small. There- 
fore, the S — state of a dot with even number of elec- 
trons is less robust than the corresponding ground state 
of a real atom. Because of the exchange interactions be- 
tween electrons residing on the last doubly occupied level, 
one of the electrons is promoted to the next (empty) or- 
bital state, and thus a triplet spin state of the dot is 
formed. Apparently, dots studied in and were in 
a triplet state in the absence of magnetic field. On the 
other hand, experiments with carbon nanotubes |^] indi- 
cate that in the absence of magnetic field the state with 
even number of electrons is a singlet. 

Application of a magnetic field results in a spin tran- 
sition in the ground state of a quantum dot or carbon 
nanotubc |^ . In the former case the Zeeman effect is neg- 
ligible because of the small g-factor in GaAs; magnetic 
field affects mostly the orbital states, making energy of 
the triplet state higher than that of the singlet. In the 
latter case, the Zeeman contribution dominates and can 
exceed the level spacing, i.e., the magnetic field induces 
the transition to a triplet state. Electron tunneling be- 
tween the dot and leads results in mixing of the com- 
ponents of the ground state. The mixing involves spin 
as well as orbital degrees of freedom and yields an en- 
hancement of the low temperature conductance through 
the dot at the transition point. This enhancement can 
be viewed as a magnetic field-induced Kondo effect JlOf . 

In Section |ll| a simple model capable of describing the 
singlet-triplet transition in Coulomb blockade systems is 
introduced. In Section III we derive an effective low en- 



ergy Hamiltonian and establish the Fermi-liquid nature of 
its ground state. The Hamiltonian is analyzed by means 



of renormalization group in Sections [V VI. Generaliza- 
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tions of the model are discussed in Appendices |C[]_ 

Section ^ is relevant for the experiments |^] with verti- 
cal GaAs quantum dots. In this case the Zeeman energy 
is very small due to an anomalously small electronic g- 
factor in GaAs. The low energy physics in this limit is 
described by a special version of a two-impurity Kondo 
model . We show that at zero temperature the linear 
DC conductance is a monotonous function of a magnetic 
field, reaching the unitary limit ~ Ae^/h at the triplet 
side of the singlet-triplet transition, and decreasing log- 
arithmically at the singlet side. The Kondo temperature 
on the triplet side of the transition falls off rapidly from 
To at the transition point down to yj^'P'''* far away from 
the transition 

The relations between Tq, rj^'P'°*^ and the Kondo tem- 
perature T'^'^ in the adjacent Coulomb blockade val- 
leys depend on the tunneling amplitudes between the 
dot and the leads and also on the intradot interac- 
tions. If the tunneling preserves the orbital symmetry, 
as it likely occurs in vertical quantum dots then 
y^ripiet ^ yodd < Accordingly, at a finite tcmper- 



ature T, 



triplet 
K 



< T < To the Kondo effect in N = 



valley of Coulomb blockade can be observed only when 
the system is tuned to the singlet-triplet transition point 

0. 

If the tunneling mixes the orbitals, the difference be- 
tween 7^j^'P'°' and T^'^ becomes less pronounced and 



K 



vanishes if the mixing is strong, see Appendix The 
relation between Tq and Tf^^ is almost independent on 
the degree of mixing but depends strongly on the intradot 
interactions, see Appendices 

In Section VI we discuss the Zeeman splitting-driven 
transition ]15|, relevant for the transport experiments 
with carbon nanotubcs . In this system, the g-factor is 
close to its free-electron value, while the orbital effects are 
suppressed due to a small radius of the nanotube. The 
low energy physics in this case is described adequately 
by the single-channel anisotropic Kondo model. Tuning 
away from the transition points is equivalent to applying 
a finite magnetic field to the conventional Kondo impu- 
rity. The DC conductance has a peak at the transition 
point at all values of T, with the corresponding Kondo 
temperature being of the same order as T^"^. 

Finally, in Section VIl we demonstrate that in a generic 



case the positions of the peaks of non-linear conductance 
are not symmetric with respect to the change of the bias 
polarity. This effect is weak and can be best seen in 
double-dot systems. 



II. THE MODEL 

We will consider a confined electron system which does 
not have special symmetries, and therefore the single 
particle levels in it are non-degenerate. In addition, we 



assume the electron-electron interaction to be relatively 
weak (the gas parameter Vg ^ 1)- Therefore, discussing 
the ground state, we concentrate on the transitions which 
involve only the lowest-spin states. In the case of even 
number of electrons N on the dot, these are states with 
5 = and S ~ \. At a sufficiently large level spacing 
5 = e+i — e_i between the last occupied (—1) and the 
first empty orbital level (-1-1), the ground state is a singlet 
aX B ~ 0. Finite magnetic field affects the orbital ener- 
gies; if it reduces the difference between the energies of 
the said orbital levels, a transition to a state with S = 1 
may occur, as illustrated at Fig. ^ Such a transition in- 
volves rearrangement of two electrons between the levels 
n = ±1. Out of the six states involved, three belong to a 
triplet = 1, and three others are singlets (S* = 0). The 
degeneracy of the triplet states is removed only by Zee- 
man energy. The three singlet states, in general, are not 
degenerate with each other. To describe qualitatively the 
transition between a singlet and the triplet in the ground 
state, it is sufficient to consider the following Hamilto- 
nian: 

Hdot = ^ndi^dns - EsS^ - EzS' +Ec{N~ Mf . 

ns 

(1) 

Here, N = „ dl^s^ns is the total number of electrons 
occupying the levels n = ±1, operator 

nss' 

is the corresponding total spin of the dot (cr are the Pauli 
matrices), and the parameters Es, Ec, and 



Ez — gdl^sB, 



(2) 



are the exchange, charging, and Zeeman energies respec- 
tively |l^. We restrict our attention to the very middle 
of a Coulomb blockade valley with an even number of 
electrons in the dot (that is modelled by setting the di- 
mensionless gate voltage M to N — 2). We assume that 
the level spacing 5 is tunable, e.g., by means of a mag- 
netic field B: 5 = 5(B). 

The lowest energy singlet state and the three compo- 
nents of the competing triplet state can be labeled as 
\S, S^) in terms of the total spin S and its z-projection 

5^ 



\l,l)=dl^/_^^\Q), 

|i,-i> = <urfLulo), 



(3) 



|o,o) = dLitrfLulo), 



<udLiJ|o), 



where |0) is the state with the two levels empty. From 
(pi), the energies of these states satisfy 
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KoS-EzS^ Ko = S-2Es (4) 



If, for example, Kq > at S = 0, the ground state of 
the dot in the absence of the magnetic field is a singlet 
|0,0). Finite field shifts the singlet and triplet states due 
to the orbital effect, and also leads to Zeeman splitting 
of the components of the triplet. As B is varied, the level 
crossings occur (see Fig. |l|). The first such crossing takes 
place a.t B — B* , for which 



A = i;|o,o> =0. 



(5) 



At this point the two states, |1, 1) and |0, 0), form a dou- 
bly degenerate ground state of the dot. 



CO 

(N 



to 




S* 5 

FIG. 1. Typical picture of the singlet-triplet transition in 
the ground state of a quantum dot. 

If leads are attached to the dot, the dot-lead tunneling 
results in the hybridization of the degenerate (singlet and 
triplet) states. The characteristic energy scale Tq associ- 
ated with the hybridization can be in different relations 
with the Zeeman splitting at field B ^ B* . 

If Ez{B*) ^ To, then the Zeeman splitting can be ne- 
glected, and at the B = B* point all four states (||) can 
be considered as degenerate. This limit adequately de- 
scribes a quantum dot formed in a two-dimensional elec- 
tron gas (2DEG) at the GaAs-AlGaAs interface, subject 
to a magnetic field @Jl^. Indeed, the orbital effect of a 
magnetic field is very strong in GaAs due to a smallness 
of the effective mass, so that a relatively weak magnetic 
field B* suffices to induce the spin transition. In addi- 
tion, the electronic g-factor in GaAs is also small. Both 
these factors contribute to the smallness of the Zeeman 
energy Ez{B*) = gdfJ'sB* compared to Tq. This case is 
studied in details in Section 

The opposite limit Ez{B*) ^ Tq is realized in car- 
bon nanotubes In nanotubes, on the one hand, the 
effect of magnetic field on orbital motion is very weak 
due to their small diameter. On the other hand, the g- 
factor is close to its free electron value g — 2, yielding an 
appreciable Zeeman splitting even in a magnetic field of 
moderate strength. This limit of the theory is addressed 
m Section 0. 

In order to study the transport problem, we need to in- 
troduce into the model the Hamiltonian of the leads and 



a term that describes the tunneling. We choose them in 
the following form: 



Hi — ^ksC^anks'^anksi S.ks - £.k - -j^gci^sB (6) 

aiiks 

Ht^ J2 iann'cl^ksdn's+il-C. (7) 

ann' ks 

Here a = R,L for the right/left lead, and n = ±1 for the 
two orbitals participating in the singlet-triplet transition; 
k labels states of the continuum spectrum in the leads, 
and s = ±1 is the spin index. The single-particle ener- 
gies ^ks include the direct effect of the external magnetic 
field on the electrons in the leads. 

In writing (|^)-(|7[), we had in mind the vertical dot de- 
vice , where the potential creating lateral confinement 
of electrons most probably does not vary much over the 
thickness of the dot . Therefore we have assumed that 
the electron orbital motion perpendicular to the axis of 
the device can be characterized by the same quantum 
number n inside the dot and in the leads. Presence of 
two orbital channels n = ±1 is important for the expla- 
nation of experiments 0, in which the orbital effect of 
the magnetic field dominates. In the case of large Zee- 
man splitting, the problem effectively becomes a single- 
channel one pSfl , as we will see in Section VI 



III. LOW ENERGY HAMILTONIAN 

We will demonstrate the derivation of the effective 
low energy Hamiltonian under a simplifying assumption 



(8) 



It has been shown previously that, despite its sim- 
plicity, Eq. (H) still allows to understand, at least quali- 
tatively, the most interesting results of the experiments 
[0,^. Some of the features not captured by (||) will be 
discussed later on. 

The advantage of (^) is that it allows us to perform a 
rotation in the R-L space pTl 



ks 



/ tR 




f Clinks \ 






\ CLnks J 



after which the cj) field decouples from the dot: 



(9) 



(10) 



nks 



The Hamiltonian of the system then acquires a form 

n = H^{4>) + H{^), (11) 

where Ha{(j)) = Y^nks ^ks4>nks'Pnks with S^ks given by 
Eq. (^) is the free-electron Hamiltonian of the 0-particles, 
and 
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(12) 



nks 



describes the '(/'-particles interacting with the dot via the 
tunnehng Hamihonian (|ic|). 

The physical quantity we are interested in is the lin- 
ear conductance of the system G. The particle current 
operator is defined as 



71 ks 



(13) 



In order to take a full advantage of the decomposition 
( ini ) , we need to express j via ip and (j) of Eq. (||) . When 
the resulting expression is substituted into the Kubo for- 
mula (see Appendix ^ for details of the calculation), it 
yields 

G = ^ 50 y de{^df/de) h^i'ImT„, (e)] . (14) 



Eqs. (|l7|)-(|T^) should be understood in a sense that the 
operators in the l.h.s. obey the same algebra as those in 
the r.h.s. The states (||) are the eigenstates of the oper- 
ators S+ and (Si • S2), while S_ and T describe transi- 
tions between the singlet state and the components of the 
triplet. Some useful relations involving these operators 
are listed in Appendix 

We can now simplify Hiip) by integrating out the vir- 
tual transitions to the states with A/" ± 1 electrons on 
the dot. This is done by means of the Schrieffer-Wolff 
transformation or, equivalently, of the second order of the 
Brilluin-Wigner perturbation theory. Using Eqs. (p^)- 
(p9|), the result of this procedure is written in terms of 
the operators Si_2 as [|lT| : 

H = Hoi^j) + K{Si ■ S2) - EzSl + Hn, (20) 

71 

Hn = J(S„„ • S+) + Vnpnn{Sl ■ S2) 

-f-^ [(s_„,„ • S_) + 2zn(s_„,„ • T)] . (21) 



Here Tns are the diagonal elements of the T-matrix for 
the Hamiltonian H{ip), defined in a usual fashion (see 
Appendix^), / is the Fermi function, ly is the density of 
states in the leads, and 



50 



2TTh \t\+t\ 



(15) 



In order to bring H{'>p) to the form convenient for fur- 
ther analysis, we explore the existence of a one-to-one 
correspondence between the low energy states of the dot 
(^ and the states of two fictitious 1/2-spins pT[ |: 



|1,-1> 
|1,0)^ 

|0,0)^ 



I T1T2), 

=^ I -1112), 
^(1 Tli2) + 



1T2)) 



(16) 



^(1 Tii2)-UiT2)) 



This correspondence allows one to represent any operator 
acting on the states (^) in terms of the spin-1/2 opera- 
tors Si ^2- By comparing directly matrix elements, one 
can check the validity of the following relations: 



Here we introduced the particle and spin densities in the 
continuum: 



Pnn' — V'i/csV'n'fc's, S„„/ — V'lfcs 2 



VY^di^dn'sV 



nSn 



2 ^ 2V2 
[(Si-S2)-l/4 + n] 



S_ + 2mT) (17) 
(18) 



Here V = J2s \S, S^){S, S'^] is a projector onto the 
low energy mul'tiplet (^, cr — ((t"^, ct^, cr^) are the Pauh 
matrices, and 



S± = Si±S2, T = SixS2. 



(19) 



(22) 



The bare values of the coupling constants are 



J = I = 2V = 2{tl+t%)/Ec. (23) 



We did not include into (p0|)-(|2l|) some terms, that 
are irrelevant for the low energy renormalization. The 
contribution of these terms to the conductance is fea- 
tureless at the energy scale of the order of Tq (see the 
next section), where the Kondo resonance develops. The 
Schrieffer-Wolff transformation also produces a small cor- 

-i?|o,o) between the 



rection to the energy gap K = E^i^gy 
states |1,0) and |0, 0), so that K ^ Kq. However, this 
difference is not important, since it only affects the value 
of the control parameter at which the singlet-triplet tran- 
sition occurs, but not the nature of the transition. 

The effective Hamiltonian (pO[)-(^l|) can be simplified 
further by excluding the magnetic field acting on the elec- 
trons in the leads. Indeed, according to Eq. (|lj), the 
linear conductance is determined by the scattering prop- 
erties of the i/i-particles in the vicinity of the Fermi level. 
It is therefore convenient to redefine the single particle 
energies in such a way that being measured from the cor- 
responding Fermi levels they are independent of the spin 
direction S,ks ^k, 



Ho 



Y^ki'iksi'nks, 
nks 



(24) 
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cf. Eq. (|12|). Most of the operators in the r.h.s. of (|2 
are invariant with respect to the redefinition of energies 
The only correction that appears is 



(25) 



where {s^„)o is ground state average of the spin density 
in the presence of the magnetic field (Note that J here 
is the bare value of the exchange amplitude). This cor- 
rection can be readily absorbed into Ez and can be cast 
in the form of a small (since i^J <C 1) tunneling- induced 
correction to the g-factor of the dot's electrons [see 
Eq. (^)], indistinguishable from other types of corrections 
Thus, not quite unexpectedly the direct Zee- 
man coupling of the magnetic field to the electrons in the 
leads has little infiuence on the properties of the system. 

There were several simplifying assumptions involved 
in the derivation of the effective Hamiltonian (pc|)-(21 ). 
Nevertheless, we believe it provides an adequate descrip- 
tion of the low temperature physics in the vicinity of the 
singlet-triplet transition. The dot interacts via tunneling 
with four species, or channels, of conduction electrons. 
At the same time, the dot behaves like a composite two- 
spin impurity. Obviously, it takes only two channels to 
screen such an impurity. By analogy with the multi- 
channel Kondo model one should expect that in the 
absence of special symmetry no more than two channels 
are coupled to the dot in the effective low energy Hamil- 
tonian. 

In the above derivation of the effective two-channel 
Hamiltonian we entirely ignored the multilevel structure 
of the dot and also adopted a simplified description of 
intradot Coulomb interaction. The presence of many 
energy levels in the dot is important at the energies 
above the single particle level spacing S, while the Kondo 
physics emerges at the energy scale well below 5. These 
levels result merely in a renormalization of the param- 
eters of the effective low energy Hamiltonian |20 . One 



only needs to consider this renormalization for deriving 
the relation between the parameters tann' of the low en- 
ergy Hamiltonian (|l]) , (||) and (^ to the "bare" constants 
of the model defined in a wide (~ ep) band. On the 
other hand, using the effective low energy Hamiltonian, 
one can calculate, in principle, the observable quantities 
such as conductance G{T) and other susceptibilities of 
the system at low temperatures {T ^ S), and establish 
the relations between them, which is our main goal. 

The two-parameter {Ec,Es) description of the in- 
tradot Coulomb interaction, see Eq. (^, is justified as 
long as Random Matrix Theory (RMT) is applicable . 
Small dots JzpUl^, containing only a few electrons, can 
not be described in the framework of RMT. Characteriza- 
tion of the interaction in such systems requires introduc- 
tion of additional parameters. Respectively, the part of 
the effective Hamiltonian which describes the interaction 



of the dot with itinerant electrons contains some addi- 
tional [compared to Eq. (^l|)] constants, see Eq. (|Cl]). 
Similar complications arise if the restriction (|^) on the 
tunneling amplitudes is lifted, see Eq. (|D2|). Although 
the Hamiltonians (D2) and (CI) appear to be different 
and more cumbersome than Eq. (|2l|), all three models 
turn out to be identical as far as one is interested in the 
low-temperature properties of the system. We will see in 
the next section that the three constants, J, V, and / 
of the exchange Hamiltonian ( |2l| ) grow in the course of 
renormalization. At a sufficiently advanced stage of the 
renormalization group procedure, these constants signif- 
icantly exceed their starting values. It turns out, see 
Appen dice s |c| a nd P , that at such a stage the Hamilto- 
nians ([D2|) and (|Cl| ) generate the same renormalization 
group (RG) flow as Eq. (|2l| ) does. 

Note that the expression of the current operator (|1J 
in terms of the fields ip Etnd (f> depends strongly on the 
manner in which tp, (f> are composed from the original con- 
duction electrons representing the left and the right leads. 
Therefore, lifting the restriction (^) results in modifica- 
tion of Eqs. (p^), which express the DC conductance G 
through the properties of V'-particles. In Section VII be- 
low we study in details the corresponding generalization 
of Eq. (|lj) onto the case of large Zeeman energy. 

Also note that Eqs. (20)-(pTl) are valid only in the 
vicinity of the singlet-triplet transition. Far away from 
the transition, the system's properties can be much more 
complicated compared to the predictions based on (pO|)- 
( pT| ) . An example of such non- universal behavior is dis- 
cussed in Appendix pj^. 

The Hamiltonian ^2^)-(|2^) is still very complex and its 
exact treatment is very difficult, if possible at all. For- 
tunately, there exists a strong argument in favor of the 
Fermi- liquid nature of the ground state of (|2^)-(21). In- 
deed, formally, (|2^)-(|2^) resembles the effective Hamilto- 
nian of the two- impurity Kondo model (2IKM) , for which 
Hn is replaced by [E2| 



2IKM 



J (s„„ • S+) + / (s_ 



(26) 



and the parameter K characterizes the strength of the 
RKKY interaction. It is known [^ that 2IKM may un- 
dergo a phase transition at some special value of K. At 
this point, the system may exhibit non- Fermi liquid prop- 
erties. However that this can happen [^^ only if H is in- 
variant with respect to the particle-hole transformation 



The two extra terms in ( |2l| ) as compared to (^6|) violate 
this invariance. Therefore, the symmetry that warrants 
the existence of the non- Fermi liquid state is absent in our 
problem. The logarithmic grow of the two extra terms 
in (j2^) at low energies []ll| precludes the possibility that 
the RG flow passes anywhere near the 2IKM non-Fermi 
liquid fixed point. 
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IV. SCALING ANALYSIS 

To gain an insight into the low-energy properties of 
our model, we apply the "poor man's" scaling technique 
p3[ . The procedure consists of a successive integration 
out of the high-energy degrees of freedom and yields a 
set of scaling equations 



dJ/dC = iy{J^ +1^), 
dl/dC = 2iyI{J + V), 
dV/dC = 2vl'^ 



(27) 



with the initial conditions (|23|). Here C = ln{S/D), and v 
is the density of states in the leads; the initial value of the 
high energy cutoff D is D — S, see the discussion after 
Eq. (^). The scahng procedure also generates correc- 
tions to K. In the following we absorb these corrections 
in the re-defined value of K. Equations ( p7| ) are valid in 
the perturbative regime and as long as 

£>> \K\,Ez,T. 

At certain value of the scale, C — Cq, the inverse coupling 
constants simultaneously reach zero: 

1/J (Co) ^ l/I (Co) = l/V (Co) ^ 0. 

This can be used to define a characteristic energy scale 
of the problem Tq through the equation 

Hd/To)=Co. 

The value of Tq is obviously non-universal. It is con- 
venient to parametrize Tq by a value of constant c in 
the expression resembling the usual definition of Kondo 
temperature: 



To = (5exp[— c/i^J]; 



(28) 



here J — J{C — 0) is given by Eq. (^31). For the special 
choice (^ of the tunneling amplitudes it was found nu- 
merically [0 that c = 0.36. It is instructive to compare 
To with the Kondo temperature j^j^'P'®*^ corresponding to 
K — a.t the triplet side of the singlet-triplet transition 
(see the next Section), 



r^"pi°' = 5exp[-l/i.J], 



(29) 



and with the Kondo temperature TJ^'*^ in the adjacent 
Coulomb blockade valleys with N — odd. In the lat- 
ter case, only one of the two orbitals rt = ±1 in the 
dot is involved in the corresponding effective Hamilto- 
nian, which takes the form of a single-channel S* = 1/2 
Kondo model. The corresponding exchange amplitude is 
Jodd ^ 4,{t\ + tj^)/Ec = 2 J, which yields 



[note that (gg) and (|3 
(i)l. In the limit lyj 



are also valid only for the model 
Eqs. (pl-d^) imply that 



To > T'^'^'^ > T, 



ntriplct 
K 



This inequality explains why the Kondo effect at the 
singlet-triplet transition point in the TV = even valley 
appears more pronounced than that in the TV = odd val- 
leys , and why the experiments with vertical quantum 
dots failed to detect the Kondo effect away from the tran- 
sition when the dot definitely was in the triplet state . 

The solution of the RG equations (^^ can be expanded 
near C ^ Cq. To the first order in Cq — C = In D /Tq, we 
obtain 



1 



A- 1 
2i^V 



(A + l)ln(i?/To), 



where 



A = 2 + \/5 w 4.2. 



(31) 



(32) 



It should be emphasized that the constant A is deter- 
mined by the properties of the solutions of (|2^) at large 
C and is universal in the sense that it's value does not 
depend on the bare values of the coupling constants. 

The solution ( |3l| ) can be used to calculate the differ- 
ential conductance at high temperature T ^ |isr| , Tg- I^i 
this regime, the coupling constants are still small, and 
the conductance can be calculated perturbatively from 
the Hamiltonian (pO[)-(^l|) with renormalized parameters 
and formula (14). This results in 



G/Go 



where Go = 4(7o and 



A 



[HT/To)r 



(33) 



A= (3^V8)(A + l)-2 1 + A+(A-1) 



0.9 



is a numerical constant. At low temperature, the RG flow 
( p7| ) terminates at either D \K\ (which corresponds to 
the relatively strong orbital effect of the magnetic field), 
or at D ^ Ez (weak orbital effect). We consider these 
two cases separately. 



V. SINGLET-TRIPLET TRANSITION 

In this section, we assume that the Zeeman energy is 
negligibly small compared to all other energy scales. At 
high temperature T ^ \K\, the conductance is given by 
Eq. (^3|). At low temperature T ^ \K\ and away from 
the singlet-triplet degeneracy point, li^l > Tq, the RG 
flow yielding Eq. ( |33| ) terminates at energy D ~ \K\. On 
the triplet side of the transition. 



Tiodd 
K 



(5exp[-l/2j/J]. 



(30) 



K « -To, 
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the two spins Si, 2 are locked together into a triplet. 
The system is therefore described by the effective two- 
channel Kondo model with S — 1 impurity, obtained 
from Eqs. (pO|)-(^ by projecting out the singlet state 
and dropping the no longer relevant potential scattering 
term: 



G/Go = F{T/Tk) 



(39) 



-fftiiplot = Hq + J ^(s„„ • S). 



(34) 



Here J is given by the solution J(£) of Eq. (p7|), taken 
at £ = £* = ln((5/|_fi'|), which corresponds to D — \K\. 
With further decrease of D the renormalization of the 
exchange amplitude J is governed by the standard RG 
equation |23 



dJ/dC = iyJ^, C>£* 



(35) 



The solution of Eq. (|35|), 

l/i^J{C)-l/i^J{£*) = C- C*, 

can be also written as 1/i/J = \ti{D /Tk) in terms of the 
running bandwidth D and the Kondo temperature 

Tk{K) = \K\ exp [-l/yJ{C*)\ . (36) 

Substituting here the asymptote of J{C*) from Eq. ( |3l| ) 
we obtain 



Tk 
To 



\K\ 



(37) 



This result is universal and is valid in the vicinity of 
the transition point. To clarify what does the "vicin- 
ity" mean, we again rely on the model (|8|), for which 
Eqs. ( p7| ) are valid for all \K\ ^ 5. Expansion of the 
solution of Eqs. (^^ near the weak coupling fixed point 
C — and the substitution of J{C*) into (36) results in 
the large- 1 if I asymptote of Tk{K): 



Tk{K) 

/jntriplct 



5 

W\ 



(38) 



is given by Eq. (||). Comparison of f^) 
|) shows that the two asymptotes match at 



\K\lTo^{5/Toy\ fi = 



l/c-2 
A- 1 



0.24, 



which gives the upper limit of the validity of Eq. (37): 
\K\/To < {S/ToY. From Eqs. (|3),(|38|) follows that the 
larger is \K\ the smaller Tk is. In fact, the universal part 
Tk{K), Eq. (|2g|), suffices to describe a fall of Tk by an 
order of magnitude for realistic values of the parameters 
for a vertical quantum dot device [ p^ . 

For a given \K\ > T,To the linear conductance can be 
cast into the scaling form. 



where F (x) is a smooth function that interpolates be- 
tween i^(a;>l) = (ttVS) (lna:)"2 ^ (q) = 1. It 
coincides with the scaled resistivity 

F{T/Tk) = p{T/Tk)/p{0) 

for the symmetric two channel 5=1 Kondo model. 
The conductance at T = (the unitary limit value) is 
Go — 4.go, see Eq. (p^). Equation ( |3^ ) remains valid not 
too far away from the transition point, where the uni- 
versal description (|20|)-(|2l|) is apphcable. At large \K\ 
the conductance exhibit additional features (an example 
is discussed in Appendix ^). However, smallness of Tk 
renders these features difficult to observe. 



-To ^ To T K 

FIG. 2. Linear conductance near a singlet-triplet transi- 
tion. At high temperature G exhibits a peak near the transi- 
tion point. At low temperature G reaches the unitary limit at 
the triplet side of the transition, and decreases monotonously 
at the singlet side. The two asymptotes merge at S> T, To- 

On the singlet side of the transition, 
K » To, 

the scaling terminates at D ^ K, and the low-energy 
effective Hamiltonian is 

-ffsinglet — Flo — -^Vy^^npnn, 

n 

where V is the renormalized scattering amplitude V{C) 
taken at C — £* — \n{S/\K\). The temperature depen- 
dence of the conductance saturates at T ^ if, reaching 
the value 



G/Go — 



B 



[HK/To)r 



B = 



37r A - 1 

TaTT 



0.5. (40) 



This result is universal in the vicinity of the transition. 

At T = Eqs. (11) and ^ predict different de- 
pendence on the parameter K which is used for tuning 
through the transition. At positive K, conductance de- 
creases with the increase of while at K ^ —To con- 
ductance G is independent of see Fig. I While the 
conductance G is a monotonous function of _ftr at T = 0, 
it develops a peak at the singlet-triplet transition point 
if = at high temperature T ^ Tq. 
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VI. TRANSITION DRIVEN BY ZEEMAN 
SPLITTING 



If the Zeeman energy is large, the RG flow (27) ter- 
minates a.t D ^ Ez- The efl'ective Hamiltonian, vahd at 
the energies D < Ez is obtained by projecting (p0|)-(|2l|) 
onto the states |1,1) and |0, 0). These states diff'er by 
a flip of a spin of a single electron (see Fig. ||), and are 
the counterparts of the spin-up and spin-down states of 
S — 1/2 impurity in the conventional Kondo problem. It 
is therefore convenient to switch to the notations |l5| 



|i,i) = |T), |o,o) = U), 



(41) 



and to describe the transitions between the two states in 
terms of the (pseudo)spin operator 



ss' 

built from the states 



ii.i>=it> 



t 

0,0) = U> 



FIG. 3. The ground state doublet in case of a large Zee- 
man splitting. The states |1, 1) and |0,0) differ by flipping a 
spin of a single electron (marked by circles). 

Among the various operators in (|20|)-(|2l|) that act on 
the states of the dot only the following have non-zero 
matrix elements between the states (p]): 

V'SIV = S' + 1/2, 7"(Si • S2)V' ^S'- 1/4, 

V2 

where V' = J2s \ ^){^\- Using these relations, we obtain 
from (|2C|)-(|2T1) 



nks 

+ E (S^ + 1/2) + VnPnn (^^ - 1/4)] (42) 

n 

-/(s+_i5-+H.c.) , 

where A = Ez — K, see Eq. (^. It is now convenient 
to transform (|4|) to a form which is diagonal in the or- 
bital indexes n. This is achieved by relabeling the fields 
according to 



(43) 



which yields 

H = Ho- A5" 



(44) 



7' „z qz 



where Hq is a free-particle Hamiltonian for a, b electrons, 
and Sa = J2kk'ss' "-is i^ss'/2) Uk's' is the spin density for 
a-electrons (with a similar definition for s^). The cou- 
pling constants in (Eih, 



Va = {J-V)/2, J, 

Vb^{J + V)/2, J'^ 



J 

J - 



2V, Ji 
2V 



2L 



(45) 



are expressed through the solutions of the RG equations 
d^) taken at £ = £* = lii{6/Ez). 

Eq. ds^ ) contains unusual terms Visf, which have a 
meaning of a magnetic field acting locally on the con- 
duction electrons at the impurity site. The appearance 
of these terms is a price of the finiteness of the mag- 
netic field, which breaks spin-rotational invariance. For- 
tunately, these terms do not change significantly prop- 
erties of the model. Their main effect is to produce a 
small correction to A, through creating a non-zero ex- 
pectation value (sf). This correction is not important, 
since it merely shifts the degeneracy point. In addition, 
ViS^ lead to insignificant corrections to the densities of 
states ||l5[. In the following, we will discard these terms. 

The scaling equations for the coupling constants in ( |4^ ) 
take the standard form [Esl 



dJz/dC = vJ\, dJ^/dC = i'JzJ±, C> C*\ 



(46) 



J'z is not renormalizcd in this order of perturbation the- 
ory, but acquires a negative correction in the next order: 
dJ'z/dC — J\J'z- In addition, initially is smaller 
than Jz (for Ez = S, for example, simply vanishes). 
Therefore the 6-dependent part of (Q) can be ignored, 
and the resulting effective Hamiltonian is that of a single- 
channel = 1/2 anisotropic Kondo model ||l5|] . 

The Kondo temperature Tr- depends on Ez- For the 
model (H), as follows from Eqs. ( ^3|) and (|4^), the ex- 
change at Ez = S is isotropic: Jz — J± = 2J (where 
J is given by (^3|)). Therefore the Kondo temperature 
corresponding to Ez = S satisfies 



T^Zccman 



Zccman 
K 



see Eq. (^0|). At lower Ez the Kondo temperature be- 
comes larger. Interestingly, the isotropy of the exchange 
in (|4j) is preserved in the regime when Eqs. (27) can be 
linearized near C — Q. The corresponding asymptote of 
Tk{Ez) is 



Tk{Ez) 

7^ Zeeman 



_5_ 

Ez 



1/2 



(47) 
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When Ez approaches Tq one can use Eqs. (31|) to cal- 
culate initial conditions for (^6|). From (^) and ( p5| ) 
follows that the difference between and Jx in this 
regime is very small: JJJ_l = a/A/2 « 1.02. Therefore, 
the anisotropy of the exchange can be neglected in the 
universal regime as well, and we obtain 



scribed by a single-channel anisotropic Kondo Hamilto- 



To \EzJ 
The asymptotes and (^) match at 



(48) 



Ez/To - {d/ToY 



3-1/c 



which separates their regions of applicability. The de- 
pendence of Tk on Ez appears to be much weaker than 
the dependence of Tk on K in the previous section, and 
practically saturates in the regime (48). 



The linear conductance is given by [cf. Eq. (14)] 



G^^go J de{~df/de) [-TTi^lmT,''{e)] 



(49) 



where 7^° is the T-matrix for the a-particles with spin 
s; 7^" does not depend on s if the term VaS^ in ( ^ ) 
is neglected. The distance to the degeneracy point A 
in Eq. ( ^ ) plays the part of the Zeeman splitting of 
the impurity spin in the conventional Kondo effect. At 
A = the Kondo resonance develops: Regardless the ini- 
tial anisotropy of the exchange constants in Eq. ( ^2|) , the 
conductance in the universal regime (when T approaches 
the Kondo temperature Tk or lower) is given by 



G = GoUT/Tk), 



(50) 



where f{x) is a smooth function interpolating between 
/(O) = 1 and f{x > 1) = (37rVl6)(lna;)-2. Function 
/{T/Tk) coincides with the scaled resistivity for the sin- 
gle channel S* = 1/2 Kondo model and its detailed shape 
is known |25|. Note that in (|50| ) the zero-temperature 
conductance Go — 2go is by a factor of 2 smaller than 
that in Eqs. (||),(||),(|^). 

At finite A ^ Tk, the scaling trajectory ( ^6|) ter- 
minates at £) ^ A. Below this energy scale the spin- 
flip scattering processes are suppressed. As a result, at 
T <C A the linear conductance is given by 



G — Gq- 



7r2/16 



M^/Tk)]' 
and is independent of temperature. 



(51) 



VII. SPIN FILTERING AND ASYMMETRY OF 
THE NON-LINEAR CONDUCTANCE CAUSED 
BY ZEEMAN SPLITTING 

The low-temperature properties of the system exhibit- 
ing the Zeeman-splitting-driven Kondo effect can be de- 



n = Hoiq>) + HoW - AS' 
2 



-Us' 



s+S" + S-S' 



') , (52) 



where s = J^kk'ss' "^Is i'^ss'/'^) i'k's'- This description is 
vahd regardless the assumption (||). The operators -0fcs 
are certain linear combinations of the conduction electron 
operators from the left and right leads. When the restric- 
tion (^) is lifted, the coefficients of such linear relation 
are s-dependent [n5|: 



</>fcs 



cos -ds sin -ffg 
— sin I?, cosi9. 



CB.ks 
CLks 



(53) 



Perhaps, the simplest possible generalization of (|^) yield- 
ing 7^ 'd-i in (|3|) istann' = tan^nn' [scc Appendix ^, 
Eq. (El)], in which case tawdg = tLs/tRs- 

The s-dependence of the coefficients in ( p3| ) may lead 
to interesting effects. We consider the linear regime first. 
It is convenient to split the current operator ( p^ into 
two parts, 

j = 30 + Sj, 

which are respectively bi-linear and quadratic in the op- 
erators ip and (f). Using Eq. (|53|) we obtain 



Jo 



d 1 
(ft 2 

d 7] 
eft 2 



^ sin(279s 

ks 



rj — s SIT? "d 



(54) 
(55) 



The operator Sj makes no contribution to the DC con- 
ductance. Indeed, since ^ si^l.^ipks + commutes 
with Ti, Eq. (|5q) can be re-written as 



(56) 



Being averaged over time, Sj yields 0: 



(Sj)^ ]im — / dt{Sjit)) 

= hm JL/s^^-to)-SHto) 



because operator 5"^ is bounded, —1/2 < {S^{t)) < 1/2. 
Thus, the contribution SG from 6j to the DC conduc- 
tance vanishes identically. 

The contribution from jo is evaluated in the same way 
as in Appendix ^ with the result 
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G^Y.9oJ de{~df/de) [^7n.lm%{e)] , (57) 



where 



5o 



2Trh 



sin2(2i3,) 



and Ts is T-matrix for i/i-particles with spin s. Equa- 
tion (^ is identical to Eq. ( |57| ) , except for s-dependence 
of (7q. Therefore, the conductance at sufficiently low 
temperature is again given by Eqs. ( |50| ) and (51) with 

At T = the electrons scatter from the dot elastically 
p^ , that is, without flip of their spin. In this regime 
the inequality g^^ =/= means that the probability of 
transmission through the dot depends on the direction 
of the spin of an incoming electron. In other words, the 
system acts as a spin filter Jisf . 

The main effect for the non-linear differential conduc- 
tance dl/dVoc consists in the violation of symmetry with 
respect to the bias sign. To see this, it is sufficient to cal- 
culate the conductance in the leading (second) order of 
the perturbation theory in J±. In this order, tunneling 
of an electron through the dot is accompanied by the 
transition of the dot between the states | t) and | |), see 
Eq. (^l|). The difference of the corresponding energies, 
_E|t\ — -E'|x)i depends on the average value of the opera- 



tor through the term JzS^S^ in the Hamiltonian (52). 
This average is not zero due to a finite DC voltage bias 
Vdc, which is described by an addition of a term 



(58) 



to the Hamiltonian (52). Using (|5^) and ( p8| ) we obtain 

(S^)= I [i^Rks^^Rks) COS^ + (clks'^Lks) sin^ §^ 



ks 

V 



where rj is given by (p5|). Therefore, the energy differ- 
ence between the states ( ^ ) acquires a bias-dependent 
correction: 



i) 



6A 



-El 
•I 



IT) 



A + ^A, 



(59) 



This correction is similar to that caused by the mag- 
netic field applied to the leads, see Eq. (|25|). But there 
is an important difference: SA in Eq. (^9|) is due to 
non- equilibrium. Obviously, the contribution to the non- 
linear conductance we discuss is suppressed if leVDcj < 
— E\iy Using Eq. ( p9|) we find for the interval of 
suppression 



^ T. ^ 

< eVnc < 



l + X 



i-x 



(60) 



Similar to the conventional Kondo problem Q higher- 
order calculation yields the differential conductance 
which peaks at the ends of the interval (^0|). Therefore 
the position of peaks is not symmetric with respect to 
the change of the bias polarity Vdc —Vdc, see Fig. |^. 

The feasibility of an experimental observation of this 
effect obviously depends on the value of the parameter 
T]. The largest possible value \r]\ = 1 is reached in a 
strongly asymmetric setup. An example of such setup 
consists of level n — +1, see Fig. ^, coupled only to the 
left lead, while the level n = — 1 is coupled only to the 
right lead. Such a control over the tunneling amplitudes 
to the individual energy levels in the dot is very difficult 
to achieve with a single-dot device. But one can instead 
think of a system of two very small quantum dots con- 
nected in series | p6[ | . If the gate voltages are tuned prop- 
erly, the double-dot system is equivalent to a single dot 
with N = even. In particular, the situation is possible 
when the wave functions of the upper and lower levels 
in Fig. ^ are localized mostly in the left and in the right 
dots respectively. The degree of the "localization" and, 
therefore, the value of the parameter ry is controlled by 
the strength of the tunneling coupling between the dots. 





A 







FIG. 4. Position of the peaks of differential conduc- 
tance at (eVoc, A)-plane. The dashed lines correspond to 
eVoc = ±A. The angle i} x/2 = r/i^Jz/i < 1 can 
be either positive or negative, depending on the sign of the 
non-universal parameter rj. 



Note that similar non-linear effects occur also in the 
case of a small Zeeman splitting, Ez <C Tq, considered in 
Section At a finite temperature and when the system 
is tuned to the singlet-triplet transition point, the differ- 
ential conductance exhibits a peak at zero bias. The peak 
is split in two when the magnetic field is tuned away from 
the transition. The location of the split peaks is given 
by the energy balance condition for activation of the in- 
terorbital scattering processes: leVocI — K + SK. Here 
SK oc eVoc is a non-universal correction analogous to 
(5 A in (|59|), which appears due to a non-zero average of 
the operator npnn in Eq. ( |2l] ) in the presence of a 
finite bias. 
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CONCLUSION 

In this paper we studied the transport through a quan- 
tum dot tuned to the singlet-triplet transition point. 
Such a transition in the ground state of an isolated dot 
can be reached by applying a magnetic field to it. The 
low temperature properties of a quantum dot attached 
to conducting leads can be described by a version of a 
two-impurity Kondo model. Depending on the value of 
Zeeman energy, two distinct physical situations are real- 
ized. 

If Zeeman energy is small, the zero-temperature lin- 
ear conductance reaches the quantum limit at the triplet 
side of the transition. At finite temperature, however, the 
conductance reaches a maximum at the transition point. 
This happens due to the extra degeneracy of the ground 
state multiplet of the dot at this point, which causes an 
increase of the characteristic temperature below which 
the conductance is approaching the unitary limit. 

If Zeeman energy is large, the ground state of the dot 
is degenerate only at the transition point. In this limit 
the model can be further reduced to the single-channel 
anisotropic Kondo Hamiltonian. The conductance ex- 
hibits a peak at the point of the transition at arbitrary 
low temperature. 

The two limits of the theory are in a good agreement 
with the results of recent experiments with vertical quan- 
tum dots and carbon nanotubes. 
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APPENDIX A: SOME PROPERTIES OF THE 
OPERATORS S± AND T 



Using 



— ^ij X^-y ^af3-ySj , it is Straightforward 



to show that the operators S± = Si ±82 and T 
obey the commutation relations 



S1XS2 



(S- 



T) - (T- S_ 
[S-,(Si.S2) 

[T,(Si.S2) 



= 4z(Si-S2). 
= -2iT, 

2 



Using these relations one can check that S_ and T are 
related by means of the unitary transformation 



Since (Si • S2) commutes with 8+, application of this 
transformation to the identity 

7 

gives at once 

7 7 

APPENDIX B: DERIVATION OF EQ. (0) 

Using Eq. (0), the particle current operator j is written 



as 



^ ^ nks 



Substitution into the Kubo formula 
e2 1 

/o 

and integration by parts yield 



G=hm--/ dte-*(bXt),j(0)]) 



G = - lim go7rt^y'lmn„fcs(w). 



(Bl) 



(B2) 



nks 



where go is given by Eq. ( |l5| ) and n„fes(a;) = 
J dte'"'n„fcs(t) is the Fourier-transform of the retarded 
correlation function 

One can express n„fcs(ci;) via the retarded Green func- 
tions of the tp and cf) particles as 

Imn„fc, {^)^\ j de[f{e + u) - /(e)] (B3) 

X ImC/°^,^(e)ImC/„fcs(e -|- w), 

where /(e) is the Fermi function, and Qnksis) is the 
Fourier-transform of 

Gnks{t) = -ie{t){{i;nks{t)ALsm)- 



In writing of Eq. (B3) we took into account that the 
Green function of the (/(-particles coincides with the un- 
perturbed (without interactions) value of Qnks , 

Therefore ImGnksje) = —n6{€ — ^ks), which removes the 
integral over e in (B3). It is convenient to represent Gnks 
in the form 



2T 



KSi-S2 



Gnks (e) = g'r^ks (e) + i^)%^ks (e)t/°fe, (e) , (B4) 
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where T„ks{e) is a diagonal element of the T-matrix. 
Since the interaction in H{ip) is local, 7^fcs(e) is in fact 
independent of k: Tnks = %is- Obviously, for e = £_ks +^ 



the second term in the r.h.s. of (B4) diverges as uj at 



small ui. Combination of Eq. (B4) with Eqs. (B2) and 



(B3), and replacement of the sum over k in Eq. (B2) by 
an integral then yield Eq. (H). 



APPENDIX C: THE EFFECTIVE EXCHANGE 
HAMILTONIAN FOR THE GENERAL FORM OF 
INTRADOT INTERACTIONS 



Two electrons occupying two single particle energy lev- 
els can form three linearly independent singlet states: 

|0,0)_i=d^_i^rfLijO), 

io,o)+i = 4i^4ijo), 

|0,0)o = -^^rfli^) |0). 

These states are the eigenstates of the Hamiltonian (|l|) 
with |0, 0)_i being the singlet state of the lowest energy. 
However, (j^) is oversimplified by the assumption of the 
validity of the Random Matrix Theory ||l^ . In a more 
general model electron-electron interactions yield non- 
zero matrix elements 

„(0,0|i/<iot|0,0)„- ^ 

for n ^ n' . As the result, the lowest energy singlet |0, 0) 
is no longer |0, 0)_i, but a certain mixture of all 3 states 
|0,0)„: 

10,0)== Cn\Q,0)n. ECn = l 

n=0,±l n 

This generalization of the dot's Hamiltonian yields ad- 
ditional coupling constants in the dot-lead exchange in- 
teraction. It turns out that such an extended model is 
identical to Eq. ( ^l|) as far as the low-temperature prop- 
erties of the system are concerned. We demonstrate it on 
a specific example JT2| , p^ : 

Ti 

Co = 0, C„ = cosa + nsina), 

v2 



The case considered in |12| , |13| corresponds to a = (or 
C_i = -Ci = l/\/2). The choice of the singlet state 
in [0 and in Eq. (||) above corresponds to a = 7r/4 (or 
C_i = 1, Ci = 0). 

Of course, the transitions between the state |0,0) and 
the triplet states |1, S*^) can be described again in terms 
of 2 fictitious spin-1/2 operators 81.2- But instead of 
(|l^)-(|l|) we have now 



-Sn,-n'{S- cos a -I- 2mT sin a). 



2 2 

SS' 

1 

+ 2^ 

P ^ di^dn'sV = nS„,n' [(Si • S2) sin(2a) - 1/4 + n] 



Accordingly, the effective Hamiltonian also modifies: 
Eq. ( pO| ) remains unchanged while Eq. ( ^l| ) is replaced 
by 

Hn = J(S„„ • S+) + V'npnn{Si ■ S2) 

+ ^(s_„,„-S_) + ^2in(s_„,„-T). (CI) 

The coupling constants here are expressed through J, /, V 
of Eq. iM) as 



Is = 1V2 cos a, It ^ lV2 sin a, V = V sin 2a 



[J is the same as in (23)]. 
The scaling equations for this Hamiltonian read 

dJ/dC = u[j^ + {ll + /I) /2] 
dIsldC = 2v {IsJ + ItV) 
dlr/dC = 21/ (ItJ + IsV) 
dV'/dC = 2vIsIt 



(C2) 



(these equations are equivalent to Eqs. (17)-(20) of |p2|). 
For a 7^ the initial values of all cou pling constants differ 
from zero. System of equations (C2) allows for a solution 
in which all the coupling constants diverge at a certain 
point C — Cq. Similar to the conventional Kondo prob- 
lem, the leading divergency is J, Is^It^V oc (£o — ^)^^ ■ 
Examination of the subleading terms allows us to con- 
clude that the solutions of (C2) and ( p7| ) coincide near 
this point if one sets Is = It = I and V' = V. Indeed, it 
follows from the second and the third equations in ( |C^ ) 
that 

^ln(/_//+) = -2i.y', /± = i(/s±/T) (C3) 
Using the leading term of the asymptote of V', 
\lvV' = (2/C) {Co-C), C > 0, 



and Eq. (|C^) we find /_//+ cx {Co - C)'^ ■ Therefore 



in the vicinity of the point Co the difference between Is 
and It can be neglected, and the resulting equations for 
J, I+,V' coincide with Eqs. ^) for J, /, V. 

So far we found a solution near some singularity point 
C — Co- We should check now that such point is reached 
indeed by solutions satisfying the proper initial condi- 
tions at £ = 0. We were only able to perform this check 
by solving numerically the system (C2). Numerical so- 



lution for a set of models defined by a in the interval 
O.OItt < a < 7r/4 indeed verified the proper divergency 
of the coupling constants. We have also found that the 
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parameter c in Eq. ( p8[ ) approaches 1/2 at small a. In 
other words, in the limit a — > +0 the difference between 
To and T^'^'^ disappears. 



The convergence of the system of equations (C2) to 
Eqs. (|2^) means that the results of the main text are 
not sensitive to the details of intradot interactions. In 
particular, Eq.(37) with A given by (|3^) holds. 



We now consider the special case a = 0, in which the 
exponent A turns out to be different. We consider a slight 
generalization of the model, allowing for the dependence 
of the tunneling amplitudes on the orbital index: ta in 
Eq. (||) is replaced by tan', see also Appendix [Note 
that for a ^ this extension does not affect the value 
of A in Eq. (^2|)]. The exchange Hamiltonian for such 
model, 



(C4) 



acquires an n-dependent exchange constant. The scaling 
equations in this case, 



(C5) 



dC 

n 

have first integral 

(J+i - J_i)/2/ = tan(/j, l-^l < 7r/2. 

The value of (p is determined by the bare value of the 
exchange constants. T he fi rst integral allows one to solve 
the scaling equations ( |C5|) exactly ||l^,|l^. Near the sin- 
gularity point Cq the solutions acquire asymptotic form 

At the triplet side of the transition at energies below 
\K\ To,Ez the system is described by the Hamilto- 
nian ( p^ ) with n-dependent exchange amplitudes J„. In 
this situation one can define two separate scales, T„ , cor- 
responding to the two exchange constants, 

r„ = \K\ exp hl/z.J„(£*)] , C* = Hd/\K\). 

At these scales the two exchange constants, j^J„, be- 
come of the order of 1. Substitution here of Jn{^) from 
Eq. then yields 



Tn 

To 



7o_ 

\K\ 



I — n sin (p 



I + n sin If 

The larger of the two scales corresponds to the smallest 
of the two exponents A,, 

minA„ = 



i2|' m 

1 — sin \ip\ 
1 + sin \ip\ 



< 1. 



The Hamiltonian ( C4 ) formally coincides with that for 
the two impurity Kondo model |22|, see Eq. (|26|) above. 



If Jn do not depend of n, this model exhibits an exotic 
non- Fermi liquid behavior ]2^ . Note that the non- Fermi- 
liquid state is extremely fragile and is also destroyed 
by potential scattering, which we have neglected in our 
derivation. The requirements a = and (p = set strict 
conditions on both the intradot interaction constants and 
on the values of the tunneling amplitudes. Finding these 
conditions is beyond the scope of the present paper. 



APPENDIX D: THE EFFECTIVE EXCHANGE 
HAMILTONIAN FOR THE CASE OF ORBITAL 
MIXING 

In order to see how our results are modified in the case 
when the orbital mixing is allowed in the course of tun- 
neling, we choose the tunneling amplitudes in Eq. in 
the form 



The rotation in the orbital space. 



1 



diagonalizes Ht - 

-ff T = ^ '^OiuC^OL^j.^dns + H.C, Van = + nt'a 
an 

and modifies Eqs. (p^) , (p^) . Next, we perform a rotation 
in the R — L space, similar to 

VRn VLn \ I CRnks 



'4^nks 

0nfcs / 2 1^ 2 \ —VLn VRn I V C-Lnks 

' V Ln ' Rn ^ ^ ^ 



which yields 



Ht^^ VnlpikJns + H.C, vl = V 



Ln + ^Rm 



We now perform a Schrieffer- Wolff transformation, which 
yields Eq. ( pO| ) with iJ„ given by 

71 

-ffn — ^71 (Sn.n ' ^+) ^ -^n (Sn.n ' 

+/\/2m (s_„,„ • T) - Vp-n,7 (Si • S2) . (D2) 
The scaling equations for this model read 
dJnIdL = y {Jl + {II + /2) /2] 

dlnldL = 2v{IV + Jnin) 

dl/dC = vYl + ^^") ^^^^ 

n 

dV/dC = v^IIn 
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In a generic case ta t'^ the initial values of all the 
coupling constants in (D3) are finite: 



J„ = /„ = 2vl/Ec, I = 2V = 2v+iv^i/Ec. 
For example, for 

= P (D4) 

one obtains 

J„-/„ = (l + n/3)Vo, / = 2y = (1 - /32) Jo, 
where Jo = 2{tl+tl)/Ec. 



System of equations (D3) allows for a solution in which 
all the coupling constants diverge at a certain point 
C — Co, with the leading divergency being of the form 
J„, /„, I,V (X {Co — C) ^^. The analysis of the subleading 
terms similar to that in Appendix ^ leads to a conclusion 
that in the vicinity of the point Cq the difference between 
J+i and J-i, and the difference between /„ and / can 
be neglected. The resulting equations are identical to 
Eqs. ( |27]) . A ccordingly, near the point Cq the solutions 
of Eqs. ( |D3| ) coincide with that of Eqs. (31) if Jn = J 
and /„ ~ I. Note that setting Jn = J and In = / in 



Eq. (D2) and performing a unitary transformation 



1 

71 



brings the exchange Hamiltonian (D2) to the for m (|2l| ). 
We solved Eqs. (p^) numerically for /? [see Eq. (D4)] in 



the interval < /3 < 0.9 (note that for the model defined 
by Eqs. (Dl) and (D4) all physical observables possess a 
symmetry with respect to the change /? — > 1//3). These 
calculations confirmed that the proper divergency of the 
coupling constants indeed occurs. 

It turns out that the relation between j'^'p^'^* and T^'^ 



(see Section tV) depends on the degree of orbital mix- 
ing: the stronger the orbitals are mixed the less the dif- 
ference between yj^'P'''* and T^'^ is. For example, for 



the model (^,(|dD T^'p'"* can be associated with the 
larger of the two energy scales Tn = 6 exp(— l/i/ J„), while 
Tgdd = 5exp(-l/i/Jodd) with Jodd - 2(1 + /32)Jo. This 
yields 



\n{S/T, 



triplet 
K 



In 



, 1 + /?^ 
■(1 + /3F 



The expression in the r.h.s. here is invariant with respect 
to the change (3 ^ 1/ (3 and its value varies between 2 for 
/3 = (or /3 = oo) and 1 for (3 = 1. Therefore, the tem- 
perature scales j^^'P'°' and T'^'^ coincide in the strong 
mixing limit /? — *■ 1. The experiments with vertical quan- 



tum dots m correspond to Tj 



triplet 
K 



odd 
K ' 



which can be 

explained only if the mixing is weak (/?^lorl//3<Cl). 

Note that the relation between Tq and T^"^ depends 
on (3 much weaker. For instance, for (3 = 1 the ratio 



In(VTo) ^ 1 
In (5/T°dd) 1 + 



0.6, 



while for /3 = this ratio equals approximately 0.7, see 
Eqs. (H) and (|5). Accordingly, for (5 > To the inequal- 
ity To ^ T^'^ holds independently of the strength of the 
orbital mixing. 



APPENDIX E: TWO-STAGE KONDO EFFECT 

In this Appendix, we consider an example of the non- 
universal features that the system may exhibit far away 
from the transition point. Consider the following gener- 
alization of (ph: 



(El) 



After a rotation in the left-right space 



tpnks 



tn \ —tLn tjin J \ CLnks 



: — ^Ln+ffin 



and a Schrieffer- Wolff transformation one obtains 
Eqs. (|2^)-(||3) with J and V in (|l|) replaced by J„ and 

Vn, 

Jn = 2Vn = 2tl/Ec, I = 2t+it^i/Ec. 

This is only a small modification as far as the vicinity 
of the singlet-triplet transition is concerned. Indeed, the 
scaling equations (|27|) are replaced by 



dJn/dC = iy{j^, + l'^) , 
dI/dC^i^lJ2(Jn + Vn), 

n 

dVn/dC = 2i//2 

If these equations are rewritten in terms of new variables 

J = (J+i + J_i)/2, V = {V+i + y_i)/2, 
J' - (J+1 - J-i)/2, V = {V+i - l/_i)/2, 

one finds that J, /, and V behave at large C according to 
Eq. (IH), all being proportional to {Cq — C)~^ , while J' 
diverges much slower, as J' cx {Cq — and V 

is not renormalized at all. 




FIG. 5. Staircase-like temperature dependence of the linear 
conductance at the triplet side of the singlet-triplet transition. 
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Interesting properties appear at the triplet side of the 
transition in the case when the exchange constants are 
very different, say, J+i ^ J_i. 

The exchange term in the Haniiltonian (^) is replaced 
by X^n i^nn ' S). One can define two separate Kondo 
temperatures 

Tn = (5exp(-l/j/J„), T+i > T_i 

corresponding to the two exchange constants J„. At 
T ^ T+i T_i the presence of the n — —1 channel 
can be neglected, so that the system in this regime is 
described by a single-channel 5* = 1 Kondo model. Ac- 
cordingly, the conductance in this regime is given by 



[5] 
[6] 



G = g+i 



^2/2 



(E2) 



where 5„ = ('^e'^/h')4.tljj^J-^. At T < T+i the scatter- 
ing in the n = +1 channel reaches the unitary limit. In 
this regime one half of the dot's spin is screened ||2^ . The 
remaining spin S* = 1/2 still interacts with the electron 
continuum in the channel n — ~1. This antiferromag- 
netic exchange coupling eventually results in the second 
stage of the Kondo effect. The corresponding asymptote 
of the conductance is 



'+1 +5-1 



3^716 
[ln(T/r_i)]^ 



r_i < T < r+i (E3) 



At T < r_i the dot's spin is screened entirely The 
zero-temperature conductance Go = g+i + g-i', the lead- 
ing finite temperature correction Gq — G (x (T/T-i)^ . 

Equations (E2) and (E3) describe a staircase- like tem- 
perature dependence of the conductance, see Fig. ||. The 
thermodynamic quantities are expected to exhibit a sim- 
ilar two-stage crossover |2^. 
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